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Abstract 

Recent coordinated observations of interplanetary scintillation (IPS) from 
the EISCAT, MERLIN, and STELab, and stereoscopic white-light imaging 
from the two heliospheric imagers (His) onboard the twin STEREO space- 
craft are significant to continuously track the propagation and evolution of 
solar eruptions throughout interplanetary space. In order to obtain a bet- 
ter understanding of the observational signatures in these two remote-sensing 
techniques, the magnetohydrodynamics of the macro-scale interplanetary dis- 
turbance and the radio-wave scattering of the micro-scale electron-density 
fluctuation are coupled and investigated using a newly-constructed multi- 
scale numerical model. This model is then applied to a case of an inter- 
planetary shock propagation within the ecliptic plane. The shock could be 
nearly invisible to an HI, once entering the Thomson-scattering sphere of the 
HI. The asymmetry in the optical images between the western and eastern 
His suggests the shock propagation off the Sun-Earth line. Meanwhile, an 
IPS signal, strongly dependent on the local electron density, is insensitive to 
the density cavity far downstream of the shock front. When this cavity (or 
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the shock nose) is cut through by an IPS ray-path, a single speed compo- 
nent at the flank (or the nose) of the shock can be recorded; when an IPS 
ray-path penetrates the sheath between the shock nose and this cavity, two 
speed components at the sheath and flank can be detected. Moreover, once 
a shock front touches an IPS ray-path, the derived position and speed at the 
irregularity source of this IPS signal, together with an assumption of a radial 
and constant propagation of the shock, can be used to estimate the later 
appearance of the shock front in the elongation of the HI field of view. The 
results of synthetic measurements from forward modelling are helpful in in- 
ferring the in-situ properties of coronal mass ejection from real observational 
data via an inverse approach. 

Keywords: Heliospheric Imaging, Interplanetary Scintillation, Multi-Scale 
Modelling 
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1.1. Interplanetary Space 

Interplanetary space can be considered as a transmission channel con- 
necting the Sun and the Earth, permeated with the ubiquitous magnetized 
solar wind flow from the Sun. The solar wind is inherently bimodal with fast 
solar wind from corona l holes and slow solar wind above coronal streamers 
(IMcComas et al.l . |2000| ) . As a consequence of the root of the spiral interplan- 
etary magnetic field (IMF) on the rotating Sun, the interface between the 
fast and slow solar winds at the same latitude gradually develops into a co- 
rotating interaction region (CIR) with increasing heliocentric distance. The 
ambient radial plasma flow in the co-rotating inhomogeneous background 
structure is frequently interrupted by coronal mass ejections (CMEs). A 
CME could undergo significant, nonlinear, and irreversible evolution dur- 
ing its interplanetary propagation, interacting with the ambient structured 
medium and other CMEs. An individual CME could have its outer magnetic 
shell stripped away to form a diffusive boundary layer by magnetic reconnec- 
tion (IWei et all 120061). be significantly pushed at its rear bo undary by a CIR 



(IDal Lago et all 120061 ). and even be entrained by a CIR (IRouillard et al. 



20091 ). The coupling of multiple CMEs in t he Sun-Earth system could result 



in a complex ejecta (IBurlaga et all 120021 ) or a shock-penetrated magnetic 
cloud (MC) (ILepping et al.l . 119971 ). Particularly, the compression effect ac- 
companying CMEs colliding can intensify the southward magnetic field and 
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23 is su bsequently responsible for major geomagnetic storms (IBurlaga et al. . 



19871 ). In terms of space weather, interplanetary space should be emphasized 



25 as a pivotal node in the development of the solar-terrestrial causal process. 
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1.2. Remote- Sensing Techniques 

The interplanetary space environment has been sampled frequently by 
observing interplaneta ry scintillation (IPS) of radio waves, beginning with 
the pioneering work of iHewish et al.l ( 11964T ) . IPS is essentially the intensity 
scintillation received at a terrestrial radio antenna as a result of a drifting 
interference pattern across the IPS ray-path, formed by density fluctuations 
in the solar wind. The intensity variance on the time-scale of 0.1 ~ 10 sec- 
onds for the micro-scale density irregularities, with a characteristic scale of 
tens to a few hundreds of kilometers, can be sequentially recorded by two 
telescopes if the baseline between tw o telescope s is basically parallel to the 
drifting speed of density irregularities (lBisil . l2006f ). Long baselines can resolve 
multiple solar wind streams crossing the IPS ray-path, which contribute to 
the integral intensity towards the Earth from the different radio-scattering 
layers along the IPS ray-path. The stream interface could be a shear and 
sliding layer between a fast stream at the high latitude and a slow stream 
at the ecliptic plane, or a compression layer within the CIR, as the IPS ob- 
servations reveal the velocity gradient and normal scintillation level for the 
sliding layer, and an inter mediate ve l ocity and enhanced scintillation level 
for the compression layer ( Bisi et al. . 2010h . T he presence of a CM E amid 
the background of bimodal solar-wind streams (IDorrian et al.l . 120081 ) can be 
identified from the IPS signal as (1) a noticeable negative lobe in the cross- 
correlation function (CCF), (2) a rapid variation of the solar wind speed, the 
CCF shape, and the scintillation level on the time-scale of hours or less. The 
polew ard deflection of th e ambient solar wind ahead of a CME on 13 May 
2005 (IBreen et al.l. I2008T) and the m icro-structure of two CMEs merging on 
16 May 2007 (IDorrian et all 120081 ) were reported by the IPS observations, 
using the European Incoherent SCATter radar (EISCAT) in northern Scandi- 
navia, the Multi-Element Radio-Linked Interferometer Network (MERLIN) 
radio telescopes in United Kingdom, and the Solar- Terrestrial Environment 
Laboratory (STELab) in Japan. The capability of the traditional and eco- 
nomical IPS technique to probe the inner heliosphere has been significantly 
improved by the current advances of long baselines of over 2000 km, the sin- 
gle obs erving frequency of a round 1.4 GHz, and dual- frequency observations 



of IPS flFallows et al.l . l2006h . 
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Modern heliospheric imaging and the traditional IPS techniques can rou- 
tinely cover the whole interplanetary space to fill in the relative observation 
data gap between the comprehensively monitored Sun and Earth. On the 
one hand, stereoscopic optical obs e rvatio ns using Heliospheric Imagers (His) 



( IHoward et al.l . l2008t lEyles et al.l . 120091 ) have been successfully realized as 



a major milestone since the launch of the t win So l ar Te rrestrial Relations 
Observatory (STEREO) spacecraft in 2006 ( iKaiserl . 120081 ). Occupying solar 
orbits approximately at 1 AU, the STEREO A leads the Earth in the west 
and the STEREO B lags the Earth in the east. Both STEREO A and B 
are separated from the Earth by 22.5° per year. Onboard either STEREO 
spacecraft, each HI instrument comprises two cameras of HI-1 and HI-2, 
whose optical axes lie in the ecliptic plane. Seen from a spacecraft, an an- 
gular distance between the Sun and a target is defined by an elongation. 
The elongation coverage is 4° ~ 24° for the HI-1 and 18.7° ~ 88.7° for the 
HI-2; the field of view (FOV) is 20° x 20° for the HI-1 and 70° x 70° for 
the HI-2; the time cadence is 40 minutes for the HI-1 and 2 hours for the 



HI-2 fjHoward et al.l . 120081 ; lEvles et al.l . |2009|). For HI imaging, the transient 
brightness of electron-scattered sunlight is blended with the stable bright- 
ness of dust-scattered sunlight, star light, and planet light. Subtracted from 
the background brightness, a white-light image could be discernible for the 
sunlight-illuminated interplanetary transients. One or two vantage views 
from HI imaging can not only recogn i ze th e background brightness of a CIR 
(IRouillard et al.l . 120081 ; ISheeley et al.l . 120081 ) , but also track the s peed, tra.jec 



tory, and shape of interplanetary transients flDavies et al.l . l2009t IWebb et al 



20091 ). An Earth-directed MC and its possible interaction with another MC 



or an incident shock could be expected to be directly imaged in interplane- 
tary space. The visibility of p otentially geoeffectiy e complicated structures 
such as MC-shock i ntera c tion (Xiong et al. . 2006al lbl) and MC-MC collision 



flXiong et all 120071 . 120091 : iLugaz et al.l . |2009|) in the FOV of an HI can give 



an early warning of a space weather event before its ultimate arrival at 1 
AU and hence at the Earth. On the other hand, STElab IPS data can be 
employed to generate a tomogr aphic reconstruction of the three-dimensional 
large-scale solar wind density (jJackson et al.l . 120031 ) , as applied to the con- 
straints of IPS data from the EISCAT and MERLIN flBreen et all 120081 ) . and 
the com parisons wi t h den sity images from the Solar Mass Ejection Imager 
(SMEI) flBisi et all 120081 ) . When an IPS ray-path lies within the FOV of 
an HI, the joint observations record simultaneous multi-scale manifestations 
for interplanetary dynamics, as demonstrated in the study of the merging 
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between two converging CMEs on 16 May 2007 ( Dorrian et all 120081 ). In 
addition, as the optical imaging from a single HI only records the horizontal 
and vertical motions of a CME, the extra depth information could be comple- 
mented by the line of sight (LOS) at a different perspective from another HI 
or from IPS. Thus, the coordinated remote-sensing observations using both 
HI and IPS cover a wide range of heliocentric distances and almost all helio- 
graphic latitudes at any time, and record the consistent multi-scale responses 
in white-light and radio wave scintillation to the interplanetary passage of a 
CME. 
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1.3. Observations and Models 

The interpretation of HI and IPS observational data relies significantly on 
theoretical modelling, since the models provide an insight into the underlying 
physical processes from its observable manifestation in the electromagnetic 
spectrum. For the integral signal along an IPS ray-path, the ambiguity of 
its LOS distribution is usually limited and removed by the addition of solar 
corona observ a tions via the iterative fitting of a kinematic solar wind model 
(jKlinglesmithl . 119971 ). This kinematic model assumes a purely radial solar 
wind flow with each stream travelling at a constant velocity, and projects 
the IPS ray-path onto a solar source surface using a ballistic mapping along 
the spiral IMF at an appropriate speed. By this method, the modes of 
solar wind occupying different regions of the IPS ray-path can be inferred 
from the distribution of dark (coronal hole) and bri ght (streamer) regions of 



the c orona onto which the projected ray-path falls (jKlinglesmithl . 119971 ; iBisi 



20061 ). Similarly, the derived distribution along the IPS ray-path constrained 
by the solar surface map can be further mapped outwards along the same bal- 
listic trajectory to the positions of Ulysess at high latitude, and Wind/ACE 
at Earth orbit, al lowing a direct c omparison with the corresponding in-situ 
observation data ( Bisi et al. . 2010l ). The intersecting points of ballistic tra- 
jectories mapped at different speeds provide a q ualitative i ndica tion of the 
location of the compression region within a CIR flBisi et all lioTbl ). Multiple 
observational data at different solar distances, longitudes, and latitudes, are 
assumed to be causally linked via ballistic mapping in this kinematic solar 
wind model, so the IPS data could be fitted wit h the extra freedom limited 
by other observational data (jKlinglesmithl . 119971 ) . For HI, the geometric in- 
fluence of the Thomson-s cattering sphere significantly a ffects the appearance 



of a CME in the FOV (jVourlidas and Howard! 120061 ; iHoward and Tappin 



20091 ). During the interplanetary propagation of a CME, the relative posi- 
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135 tion between the CME and the Thomson sphere changes continuously and 

136 different parts of the CME would be imaged by the HI at different times. 

137 For simultaneous stereoscopic imaging, the Thomson-scattering responses of 

138 the HI-A and HI-B could be quite asymmetrical. Moreover, the coexistence 

139 and possible interaction between the co-rotating background plasma in a 
CIR and the outwardly moving CME, can further complicate the interpre- 
tation of white-light images. Because of the intractable nature of analytical 

142 methods, the numerical model can assist by providing global context and 



140 
141 



144 



143 hints of what can and cannot be observed in HI images (jOdstrcil and Pizzo 



20091 1 . Meanwhile, the time-dependent solution from forward modelling with 

145 a numerical heliospheric model can give the profile along any IPS ray-path, 

we Thus, the local IPS data are assimilated into the global numerical model in 

147 a self-consistent way. Further, if an IPS model for the micro-scale scatter- 

148 ing process is logically coupled with a heliospheric model for the macro-scale 

149 magnetohydrodynamic (MHD) context, the synthetic observation signatures 

150 of both HI and IPS can be derived from the multi-scale model and be linked 

151 to the interplanetary dynamics. Motivated by the integration of observa- 

152 tional data as irrefutable evidence, and theoretical modelling as a convincing 

153 interpretation, we have undertaken a preliminary study in this paper of an in- 

154 cident shock propagation within the ecliptic plane using a newly-constructed 

155 multi-scale model. The propagation and evolution of an interplanetary CME 
we can be better monitored and understood by the joint efforts of an inverse ap- 

157 proach from the remote- sensing observations and a forward modelling from 

158 the numerical simulations. 

159 In this paper, the synthetic observation signatures of white-light imaging 

160 and IPS are forwardly modelled for the propagation of an interplanetary 
lei shock within the ecliptic plane. We present the multi-scale numerical model 

162 in Section 121 describe the individual observations from the His in Section [3] 

163 and IPS in Section HJ analyze the correspondence and coordination of the 
«4 simultaneous white-light data and IPS data for the same spatial position 
165 of interplanetary space in Section |5l and summarize this paper and discuss 
we the potential of using the coordinated HI and IPS observations to identify a 
167 possible longitudinal deflection of a CME/Shock in Section 



lea 2. A Multi-Scale Numerical Model 

169 The observation signatures of white-light brightness and radio-wave scin- 
no tillation are physically described by our newly constructed numerical model, 
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an integration of a numerical MHD model for the macro-scale driver and 
an IPS model for the micro-scale response. This multi-scale model is im- 
plemented in the following steps to causally couple the multiple physical 
processes at dramatically different temporal and spatial scales. 

The global MHD mo del f or int erplanetary space, governed by a set of 
ideal MHD equations (c.f. lHul . ll998l ). is numerically solved by a mathematical 
shock- capturing al gorithm of Piecewise Parabolic Method with a Lagrangian 
remap (PPMLR) Jcolella and Woodward! Il984l ; iHu et al!l2007f ). The eclip- 
tic plane from 0.1 to 1 AU is discretized into a mesh of 200 x 360 grids with 
its radial spacing 0.0045 AU and longitudinal spacing 1°. The background 
solar wind is prescribed by the inner boundary conditions at 0.1 AU: number 
density of 8.2 x 10~ 4 m -3 , magnetic field of 250 nT, solar rotation speed of 
2.9 x 10 -6 rad s _1 , a temperature of 4.65 x 10 5 K, and a bulk-flow speed 
of 375 km s _1 . In this model, the background solar wind is uniform along 
the longitudinal direction. Then, the initial equilibrium of interplanetary 
space is disturbed by a fast MHD shock from solar eruption, characterized 
by the shock nose at longitude of 45°, the shock front width of 30°, the ini- 
tial shock speed 1150 km s _1 , the total pressure ratio of 18 across the front 
discontinuity, the disturbance duration of 2 hours. The introduction of this 
shock into the simulation domain is numerical ly realized by the m odification 



of the inner boundary conditions at 0.1 AU f Xiong et al. . 2006al Jbl). Thus 



the propagation and evolution of a shock wave through interplanetary space 
is quantitatively described in a macro-scale by the MHD model. 

The synthetic white-light imaging is generated fro m MHD simulation data 
by the well-established Thomson -scattering theory fjVourlidas and Howard! 
20061 ; iHoward and Tappin! 120091 ). The location of the maximum scattered 



sunlight lies on the Thomson-scattering sphere which is centred half-way be- 
tween the Sun and the observer with its diameter equal to the Sun-observer 
distance. The geometry factor of the Thomson-scatte ring sphere should be 
included to correctly interp r et the heliospheric imaging (IVourlidas and Howard 



20061 ; IHoward and Tappin! . |2009i ). Both the total electron number and its 



position relative to the Thomson-scattering sphere si gnificantly affect the 
bright ness in interplanetary space. As pointed out by IHoward and Tappin 
(120091 ). (1) though the Thomson scattering itself is minimized on the Thom- 
son surface, the scattered sunlight is maximized on the Thomson surface; (2) 
the scattered sunlight is maximized simply because it is at the point along 
any LOS that is closest to the Sun, where the incident sunlight and den- 
sity are greatest; (3) the scattered intensity becomes more spread out with 
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distance fr om the Thomson surface. So the importance of the Thomson sur- 
face in the I Vourlidas and Howard! (120061) method is somewhat de-emphasized 



in the Howard and Tappinl 



20091 ) method. However, in this paper, only 
Vourlidas and Howard! (12006!) method is adopted for simplicity. The more 



sophisticated method of iHoward and Tappinl (120091 ) will be considered in our 
model in the future. For an interplanetar y electron scattering photosp heric 



light, the geometry sketch is illustrated by I Vourlidas and Howard! (120061 . Fig- 
ure 1), the Thomson-scattering algorithm is available as an IDL procedure 
"eltheory.pro" in the solar-soft library under the SOHO/LASCO directory. 
Hence, the propagation of interplanetary disturbances in our model directly 
provides their manifestation in the synthetic white-light images. 

The IPS model for the local intensity modulation of radio w aves makes use 



of a Born approxima tion to the general weak-scattering theory (ITatarski et al. 



19931 ; Fallows! . 1200 ll ). The spatial fluctuation of the local density irregular- 



ities at a micro-scale of about 200 km is conveyed by the ambient solar 
wind flow, and consequently introduces a scintillation pattern in the (x, y) 
plane perpendicular to the IPS ray-path along the z direction. Here (x, y, 
z) is a cartesian coordinate system centred on the Earth. The amount of 
scintillation introduced by density irregularities in the solar wind varies ac- 
cording to the square of density fluctuations 5N% (jSalpeterl . 119671 ). As the 
variation of SN e with heliocentric distance is not well determined from ob- 
servations, it is a com mon practice in IPS stud i es to assume that SN e varies 
as N e , as suggested by Houminer and Hewish ( 1972I ). The intensity P&i in 
a total spectrum of the spatial wave vector (k x , k y ) is merely a linear su- 
perposition of all contributions from every thin scattering layer along the 
IPS ray-path connecting the Eart h to a remote r adio s ource, as described by 
PAi(k x ,k y ) = J P' AI (k x , k y , z)dz (jKlinglesmithl . 1 19971 ). For each scattering 
layer with its thickness dz at distance z, a linear relation between the radio 
intensity scintillation P' AI (k x ,k y , z) and the spectrum of electron density ir- 
regularities, Pn,., is given by the following set of mathematical expressions 
with the notations of the electron number density N e , the classic electron 
radius r e = 2.82 x 10~ 5 m, the observing wavelength A, the spatial wave 
number k, the Fresnel radius rf, the spectral exponent a, the dissipation- 
associated inner scale of the sharp drop of spectral power k c , the wave vector 
parallel (perpendicular) to the magnetic field k\\ (k±), and the axial ratio 
of anisotropy degree AR at a heliocentric distance r with its reference of 
AR = 8 at r = 5 solar radii ( Coles and Harmon . 19891 ; Klinglesmith . 1997 ; 
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MasseyUl998f l: 



P' AI (k Xj k yj z) = 87r(r e A) 2 sin 2 (— )P Ne (k x ,k y ,k z = 0, z)dz (1) 

r, = VTz (2) 

P Ne (k x ,k y ,k z ,z) cc P Ne (k h k±) cc N 2 e - k 2j AR2 ^ a/2 (3) 

AR = (AR - 1) • (r /r) L3 + 1 (4) 

247 The cylindrically-symmetric coordinate system (k\\,k±) with its axis along 

248 the local magnetic field is transferred from the coordinate system (x, y, z) de- 

249 fined with respect to the IPS ray-path. In equation (JTJ, the term sin 2 ( ^ h [^ z ) 

250 is a high pass Fresnel filter with its Fresnel radius rj = \/Xz. Such a radius 

251 Tf determines the maximum scale of the irregularities, at which the ampli- 

252 tude fluctuation can be received at the Earth. When the spatial spectrum of 

253 intensity PAi{k x , k y ) carried by an anti-sunward speed V is drifted across an 

254 IPS ray-path with its intersection angle (90° — 9), only the speed component 

255 perpendicular to the IPS ray-path |V drift | = |V| ■ cos# is detectable by a ter- 

256 restrial radio antenna. The geometry factor cos 8 varies along the entire IPS 

257 ray-path. Moreover, the nearly identical diffraction pattern of radio signa- 

258 tures can be sequentially received by two fixed terrestrial antennas separated 

259 by a baseline b. If the baseline b on the Earth is approximately parallel to 

260 the drifting direction V drift of micro-scale density irregularities in interplane- 
tary space, the scintillation patterns at the two telescopes will be correlated 
with some time lag r. At any scattering layer z, the spatial correlation func- 
tion between two antennas i?' 12 (b, r, z) = -R' 12 (S = b — V dr 

ift ■ 7") 2>) can be 

264 converted from the spatial spectrum P' AI (k x , k y , z) via a Fourier transform. 

265 Specifically, the spatial-to-temporal conversion is merely a cut in the spatial 

266 correlation function along the direction of b — V drift • r. When the base- 

267 line b is zero, the Cross- Correlation Function (CCF) between two antennas 

268 is degenerated to the Auto-Correlation Function (ACF) for either antenna. 

269 For a single scattering layer, the CCF is simply derived by shifting the ACF 

270 at a time lag |b|/|V drift |. Inversely, the drifting speed V drift can be inferred 

271 from the time lag r, which is the principle of IPS estimation for solar wind 

272 outflow. The drift speed in the CCF of IPS signals is the manifestation of 

273 the flow speed of the local density irregularities. As the amplitude of Alfven 

274 turbulence is much smaller in interplanetary space compared with the inner 
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solar corona, IPS speed is li kely to be close to th e bulk plasma flow speed, at 
least in the slow solar wind (jKlinglesmithl . 119971 ). In the fast so lar wind close 
to th e Sun, IPS speeds may overestimate solar wind speeds (jKlinglesmithl . 
19971 ). However, as a remote-sensing technique, an IPS signature involves 



the contributions from all of the scattering layers along its IPS ray-path. 
The ability to resolve the distribution of solar wind velocities from the cross- 
correlation of IPS measurements at two sites depends upon the differences 
in time-lags between the maxima in cross-correlation produced by different 
streams, and thus on the baseline length b and the Fresnel radius 77 (Equa- 
tion [2]). In general, l onger baselines improve the ability to resolve streams 
of different velocities (IBreen et all 120081 ; iBisi et all |2010| ). provided that the 
time interval is long enough to derive well-defined scintillation spectra, when 
the observing geometry is suitable for cross-correlation analysis. In our IPS 
model, the observing wavelength (frequency) is 21 cm (1420 MHz) and the 
baseline b is 2000 km. For the radio wave at 1420 MHz, the maximum scale 
of the irregularity from a scattering layer at a distance of z — 1 AU is 177 
km. As the local IPS signature is influenced by the global parameters of 
bulk-flow speed, electron density, and magnetic-field orientation, synthetic 
IPS data can be hierarchically generated from MHD simulation data. 



294 3. White-Light Imaging 

295 The propagation of an interplanetary shock can be continuously tracked 

296 at a macro-scale by white-light imaging of the inner heliosphere. Synthetic 

297 brightness images of interplanetary disturbances can be generated from a nu- 

298 merical MHD model using the Thomson-scattering principle. In our model, 

299 an incident shock initially launched 45° west of the Sun-Earth line, is charac- 

300 terized by a longitudinal width of 30° along its front and a speed of 1150 km 

301 s _1 at its nose. The time-series evolution of the shock is shown in Figured] for 

302 density n and Figure [2] for radial speed v r . A noticeable trailing cavity with 

303 low density is formed and expands as the shock front propagates out from the 

304 Sun. Between this so-called density cavity and the shock front lies the sheath 

305 region. Across the shock front towards the sheath, the spiral IMF lines are 

306 significantly distorted and compressed, and the bulk flow speed v r , number 

307 density n, plasma temperature are abruptly enhanced. Hence the sheath 

308 downstream of a shock front is a readily observable target in interplanetary 

309 space. White-light images are simultaneously simulated for the twin His at 

310 points "A" and "B" in Figure [H With the longitude of 45° beside the Earth 
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311 at 1 AU, any HI has the FOV from 6° to 60° in the elongation. From the 

312 combined views of HI- A and HI-B, the Sun-Earth line is completely covered, 

313 and interplanetary space is routinely monitored. However, the sensitivity 
3u of HI for a remote plasma parcel depends on not only its innate electron 

315 number, but also its heliographic position. Such a geometry dependence for 

316 white-light imaging arises from the Thomson-scattering process, the working 

317 principle of the HI instrument. The Thomson-scattering effect is strongest 

318 at the sphere marked by a dotted white circle in Figure HJ when the in- 

319 coming direction of an incident photospheric photon is perpendicular to the 

320 outgoing scattered light towards the receiver. For the shock studied in this 

321 paper, the Thomson-scattering effect is very weak for the white light received 

322 at HI-A, as the shock propagates along the diameter of Thomson sphere of 

323 HI-A. Moreover, only the western flank of the shock front is within the HI- 

324 A's FOV, and the distance between this shock flank and Thomson sphere is 

325 large due to the finite front width. In addition, because the LOS from HI-A 

326 penetrates both the sheath with high density and the trailing cavity with 

327 low density, the integral signal of remote sensing can only give an average 

328 effect. Therefore, the brightness of this shock in HI-A's FOV is very faint, as 

329 shown in Figure [31 The two-dimensional time-elongation image is assembled 

330 from a series of one-dimensional slices taken at different snapshots within 

331 the ecliptic plane. This time-elongation format, widely used in the analysis 

332 of real observational data, has been shown to be very e ffective at revealin g 



333 the evolution of solar ejecta from white-light imaging (IDavis et al.l . 120091 ). 

334 Moreover, as the background brightness 1$ abruptly decreases away from the 

335 Sun, the relative-brightness enhancement (/ — Iq)/Io is adopted to highlight 

336 the brightness deviation from the initial steady state. The relative brightness 

337 in the HI-A's FOV is far less than 0.1 until the shock front approaches and 

338 crosses the Thomson sphere near 1 AU. Hence the shock, heading towards 

339 the HI-A, is essentially invisible within the HI-A's FOV. As a dramatic con- 

340 trast, a bright diagonal streak is conspicuous in HI-B data with its relative 

341 brightness up to 0.9. This bright streak is immediately followed by a dark 

342 streak with a sharp boundary between them. The bright and dark features 

343 in Figure [3b correspond to the relative positions of the sheath and cavity 

344 in the HI-B FOV in Figure [TJ The sheath and cavity are separately imaged 

345 via different LOSs from HI-B. The shock nose and the eastern flank are con- 

346 tinuously detected by one varying narrow LOS band from the HI-B, whose 

347 elongation is 24° ~ 27° at 19 hours, 28° ~ 32° at 23 hours, and 35° ~ 41° 

348 at 33 hours. The asymmetry between the white-light images of the HI-A 
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and HI-B results from the initial deviation of shock propagation from the 
Sun-Earth line. 

With some assumptions or modelling, the spatial position of a local high- 
density structure could be inversely inferred from its contribution of the 
Thomson-scattering emission to a brightness feature in the white-light im- 
ages. Once the Thomson-scattering source is pinpointed, the total electron 
number inside the source region can be accurately calculated from the bright- 
ness corrected by the geometry factor of the Thomson-scattering sphere. 
Without an additional observation limitation, the intersection point between 
the Thomson sphere and a LOS may be simply assumed to be the source 
for the corresponding bright pixel in the white-light image. This rough as- 
sumption of an interplanetary scattering source lying on the Thomson sphere 
suffers from a significant underestimation of total electron number. The un- 
derestimation degree has been studied for the radial propaga tion of a single 
electron at various longitudes ( iVourlidas and Howardl . 12006k Figure 5). Di- 
rectly generating the behavio r of a single propagating electron to a CME, 
Vourlidas and Howardl ( 2006 ) found from their qualitative model that the 



mass underestimation (1) exceeds a factor of 2 for a limb CME at elongations 
larger than 60°, (2) also exceeds a factor of 2 for a halo CME at elongations 
smaller than 20°, and (3) is never off by more than 20% for a CME, prop- 
agating along intermediate longitudes (~ 40°), even at extreme elongations. 
As a contrast to the results from heliospheric imaging, the mass underesti- 
mation from the coronagraph for heliocentric distances of less than 30 solar 
radii is less than 50%, even when a simpler assumption is made that the elec- 
tron source is solely at the so-named plane of sky perpendicular to the Sun- 
observer line ( IVourlidas et al.l . |2000| ). Moreover, the real three-dimensional 
density distribution is more complex for an interplanetary propagating CME 
against the background of the ambient bimodal solar wind, which is difficult 
to retrieve from white-light imaging. A more realistic investigation has to 
resort to numerical simulations, particularly for the explanation of a practi- 
cal observation event. Though the outline of an interplanetary CME in the 
white-light image can be easily identified by the excess brightness of CME 
images subtracted from a pre-event background image, the conversion from 
such an excess brightness to the actual mass is not straightforward. Simul- 
taneous imaging from two vantage points of the twin STEREO spacecraft 
can significantly improve the capability of identifying the spatial location of 
electron source for an Earth-directed halo CME. However, when a front-side 
CME propagates off the Sun-Earth line and fully enters the Thomson sphere 



12 



387 of one HI, only the other HI can discern the CME-driven shock front, as 

388 demonstrated in Figures [T] and [3] of this paper. Therefore, in order to reduce 

389 the ambiguity of interpretation for the observation data of STEREO HI-A 

390 and HI-B, some additional observation techniques are necessary to provide 

391 some further observational constraints. 



392 4. Interplanetary Scintillation Signal 

393 Measurements of radio scintillation can provide a way of probing the 

394 physical processes at a micro-scale of ~ 200 km, such as the electron density, 

395 bulk flow speed, magnetic-field direction, and level of Alfvenic turbulence. 

396 Because an IPS signal is roughly proportional to the square of electron density 

397 (equation [3]) and an HI signal depends linearly on the electron density, the 

398 IPS signal is even more sensitive to high density regions of solar wind. When 

399 the tiny source of a strong IPS signal is located in the bright domain of HI 

400 imaging, the IPS signal is then confirmed to be the micro-scale manifestation 

401 of a macro-scale interplanetary transient. Thus, the HI and IPS data can be 

402 correlated and complement each other. 

403 In this model, the synthetic IPS data are generated from the distribution 
of the MHD data along the IPS ray-path of elongation 30°. For the afore- 



404 



405 mentioned incident shock, its front, sheath, and cavity subsequently cross 



406 



the IPS ray-path at 19, 23, and 33 hours respectively (Figures [T] and [2]). The 

407 response of radio scintillation to the shock passage is shown in Figure and 

408 compared with the background state in Figure H] to highlight the differences 

409 between them. For the ACF at 19 hours (Figure [5b), the negative dip occurs 

410 as a result of oscillation of the ACF beside its central peak. By contrast to 

411 the ACF for the background solar wind in Figure Hb, the angle between the 

412 IPS ray-path and the local IMF line is found to be changed from 0° to almost 

413 90°. The rotation of IMF lines in the scattering source of an IPS signal gen 



415 
416 
417 
418 
419 
420 



422 



414 erally means the passage of a CME across the IPS ray-path (IDorrian et al. 



20081 1. consistent with the global magnetic- field configuration from the MHD 



model (Figure [Tfej). In addition, the relatively small amplitude of the negative 
dip is ascribed to a small axial ratio in the micro-scale interplanetary irreg- 
ularities (equation H]), as the anisotropic distribution of density irregularities 
with respect to the magnetic-fiel d line is dramatically reduced betw e en th e 



corona and interplanetary space ([Armstrong et al.l . Il990l ; I Grail et all 119971 ) 



421 Further, the multiple streams across the IPS ray-path could be recorded as 
corresponding multiple peaks in the CCF. Given an IPS baseline parallel to 
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423 the solar wind flow direction and long enough to separate multiple peaks in 

424 the CCF, the time lag for each peak in the CCF can be read to infer the flow 

425 speed of its corresponding stream. As an instant response to the arrival of 

426 the shock front, the negative bay in the CCF is obviously intensified, and the 
42? time lag of the CCF is reduced from 5.1 seconds (Figure HE) to 2.9 seconds 

428 (Figure Eli)- With an IPS baseline of 2000 km, the bulk flow speed perpen- 

429 dicular to the IPS ray-path v t can be calculated. Across the shock front, 

430 the flow speed Vt is abruptly increased from 392 to 690 km s" 1 . The shock 

431 strongly disturbs the interplanetary medium as it passes by. At 23 hours, 

432 double peaks appear in the CCF. Their amplitudes are 0.24 at —3.3 seconds 

433 and 0.14 at —3.8 seconds in Figure [Sh, far less than the previous amplitude 

434 of a single peak of 0.37 in Figure (SU. Two distinct flows with v t = 606 km s" 1 

435 and 526 km s" 1 (Figure Eh) coexist along the IPS ray-path, which correspond 

436 to the sheath and flank of this shock (Figure [Tb, [2b, [5b, and[5f). Though the 

437 shock flank occupies a smaller section along the IPS ray-path, it has a higher 

438 density. In other words, the smaller number of scattering layers is largely 

439 offset by the stronger scintillation level in each scattering layer. The total 

440 scintillation signal from the flank is comparable to that of the sheath. As the 

441 shock moves on, the sheath is replaced by the cavity along the IPS ray-path 

442 (Figure [lb and [2b). However, the cavity with very low density contributes 

443 little to the IPS signal, and is hence ignored. Only the flank with v t = 385 

444 km s -1 could be captured by the CCF at the time lag of 5.2 seconds (Figure 

445 EJ). The feedback of IPS measurement (Figure [5]) to the HI imaging (Figure 

446 [3]) could discern from which depth of HI LOS the major brightness comes. 



447 5. Coordinated Observations of Heliospheric Imaging and Inter- 

448 planetary Scintillation 

449 The continuous increase of elongation for a bright pattern in an HI image 
is the manifestation of a shock front driven by an outwardly propagating 
CME from the Sun. The movement of the bright front is at the fast shock 

452 speed v shock . A fast shock is formed as a result of intersecting of characteristic 

453 lines of fast magnetosonic wave v r + Cf upstream and downstream of the wave 

454 front. The fast shock, behaving as a sharp discontinuity, is faster than the 

455 bulk flow speed v r , and slo wer than the fast magneto sonic wave v r + Cf 



450 
451 



456 just downstream of its front fjjeffrey and Taniutil . Il964l ). For a case of the 



457 incident shock in this paper, these characteristic speeds just downstream of 
shock front are shown in Figure [7b. The furthest point within the shock 



458 
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nose from the Sun, defined as a shock aphelion, is continuously tracked and 
presented in Figure [7^. The slope of the shock aphelion in Figure [Th, is the 
shock speed f shock in Figure [7b. Bounded by the bulk flow speed v r and fast 
magnetosonic speed v r + Cf, the shock speed v shock is gradually decreased from 
1150 km s _1 at 0.1 AU to 670 km s _1 at 1 AU during the transiting time 
of 45 hours. As demonstrated in Figure [7b, the bulk flow speed v r is far 
greater than the fast wave speed Cf because of the supersonic solar-wind flow 
and the radially- decreased IMF strength. Thus, the bulk flow speed v r just 
downstream of a shock front is a reasonable approximation to the true shock 
speed in interplanetary space, as demonstrated in this numerical case with 
the underestimation being less than 10%. Within observational accuracy, 
there would generally be a speed match between a global bright front in a 
white-light image and a local density irregularity in an IPS signal, if the 
IPS ray-path lies within the FOV of the white-light imaging. For instance, 
on 16 May 2007, two converging CMEs were merged to form a discernible 
front in the STEREO HI-A observations, whose speeds in the plane of sky 
were 325 km s _1 and 550 km s " 1 from the HI imag e, and 420 ± 10 km s _1 
and 520 ± 20 km s -1 from IPS (IDorrian et al.l . 120081 ). The speed agreement 
supports the coincidence of white-light data and IPS data for the same solar 
eruption event. 

With a simultaneous observation of IPS as an additional limit, the three- 
dimensional anti-sunward movement of a shock front could be quantitatively 
linked to its manifestation as an outwardly-moving bright front in a two- 
dimensional white-light image. Continuous white-light imaging presented in 
a time-elongation format (Figure Et) can give the receiver-associated angular 
speed Q of the bright pattern. However, the angular speed Q for each plasma 
parcel along the ray-path of each LOS is quite different. As a result, the 
plasma parcels imaged earlier by one LOS would be cut later by a series of 
adjoining LOS rays. The profile of angular speed Q along a LOS is unknown 
from the practical optical imaging, as remote sensing only gives the final 
integral effect along each LOS. As a contrast, the profile of various parameters 
along every LOS is available from a numerical model. For the numerical case 
of this paper, the radial speed from the Sun, v r , and the angular speed 
relative to the HI-B receiver, Q mB , of each plasma parcel with its relative 
brightness contribution, n ■ r| IB , are resolved along each LOS, r HIB (Figure 
[6]). As the shock is far away from the HI-B, the maximum angular speed, 
f2 HIB , for the effective brightness contribution would be shifted from the shock 
nose (Figure [6b) to the eastern shock flank (Figure [6}). The match of a shock 
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497 aphelion in interplanetary space to the outermost brightness elongation in an 

498 HI image only happens at the near-Sun distance. For this numerical case, 

499 such an elongation deviation occurs at 20 hours (Figure [7b), corresponding to 

500 the radial distance of 0.57 AU (Figure Uh)- At 19 hours, with the shock front 

501 being initially cut by an IPS ray-path (Figure [lb,), one specific component 

502 of the shock speed perpendicular to the IPS ray-path is approximated to be 

503 690 km s _1 from the IPS observations (Figure EH). Under the assumption 

504 of radial propagation, the three-dimensional shock speed is then calculated 

505 to be 711 km s" 1 by its projection measured with the IPS signal. Because 

506 the intensity of an IPS signal roughly depends on the square of electron 

507 density (equation [3]), and the background electron density drops as result 

508 of the solar wind expansion, the IPS scattering source is very close to the 

509 so-called "p" point in the literature, the closest point along the IPS ray- 

510 path to the Sun. Further, the position of the plasma parcel detected as 

511 an IPS scattering source could be calculated by the intersection between 

512 the IPS ray-path from the Earth and the FOV of the most brightness from 

513 the HI (Figure [lb,). With the derived position and speed, the imminent 

514 trajectory of the plasma parcel of the IPS scattering source at 19 hours 

515 could be predicted, whose manifestations in the radial distance from the Sun 

516 and the elongation from the HI-B are shown as black dashed lines in Figure 

517 [7b and [7b, respectively. In terms of the radial distance and the elongation, 
sis this particular plasma parcel follows the shock aphelion. Moreover, in Figure 

519 [7b, the slope of the plasma parcel is obviously smaller than that of the bright 

520 pattern. The lag of the plasma parcel in the elongation of HI-B (Figure 

521 [7b) is ascribed to the relative distance from HI-B. Located at the eastern 

522 flank of the shock front, the plasma parcel detected by an IPS signal at 

523 19 hours is further away from HI-B than other parts at the eastern flank 

524 (Figures H^, and [2b,). For this plasma parcel, the longer distance from the 

525 HI-B, r HIB , slows down the relative angular speed, f2 HIB , as demonstrated in 

526 Figure [6b-c. According to this numerical case, the predictable appearance 

527 in the elongation of an IPS-detected plasma parcel could serve as a lower 

528 limit for the outermost elongation of an outwardly propagating bright front 

529 in the heliospheric imaging. By the coupling of white-light imaging and IPS 

530 signal, the interplanetary process of a CME/shock can be better described 

531 and understood. 
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532 6. Summaries and conclusions 



533 The observational signatures of white-light imaging and IPS for the propa- 

534 gation of an interplanetary shock through the ambient slow solar wind within 

535 the ecliptic plane is analyzed via forward modelling from a newly-constructed 

536 multi-scale numerical model. This numerical model directly linking inter- 

537 planetary dynamics to observational signatures is summarized as a flow chart 

538 in Figure [H A shock front can be sharply captured and continuously tracked 

539 within the FOV of white-light imaging, once being near the surface of the 

540 Thomson sphere of a receiver. The stereoscopic imaging from two spacecraft 

541 beside the Earth can well monitor an Earth-directed Halo CME, when the 

542 FOVs from these two vantage points are simultaneously focused towards the 



544 
545 



567 



543 Sun-Earth line. As demonstrated by iDavis et al.l (120091 ) for a typical Earth 



directed CME on 13 December 2008, the CME viewed as a halo CME in the 
coronagraph image from the Earth was symmetrically imaged by the HI-A 

546 and HI-B onboard the two STEREO spacecraft, and was predicted about its 

547 speed and direction at least 24 hours before its arrival at the ACE spacecraft 

548 near 1 AU. But, if a front-side CME propagates off the Sun-Earth line, the 

549 records of two His are asymmetric. Probably, the CME is invisible to one 

550 of the His, if fully entering its Thomson sphere. In this case, the Thomson- 

551 scattering source in interplanetary space is difficult to locate on basis of the 

552 white-light imaging from the remaining HI. However, the ambiguity in locat- 

553 ing the three-dimensional spatial position from the two-dimensional bright 

554 front can be more or less relieved with the aid of additional IPS data, if the 

555 IPS signal and white-light imaging are coincident for the same CME event. 

556 Being cut by an IPS ray-path, the high density-region downstream of a shock 

557 front can be measured in terms of its bulk flow speed. When both LOSs of HI 

558 and IPS simultaneously target the shock nose, the local plasma parcel at the 

559 intersection point can be estimated about its spatial position and flow speed 

560 at that time. With the assumption of radial propagation, the plasma parcel 

561 can be predicted about its trajectory. As the bulk flow speed just down- 

562 stream of a shock front is very near to the shock speed in interplanetary 

563 space, the trajectory of the plasma parcel is a slower limit for the marching 

564 shock front. Therefore, the appearance of the predicted plasma-parcel trajec- 

565 tory in the HI FOV could serve as a lower limit for the outermost elongation 
see of an outwardly propagating bright front in the white-light imaging. 

As the most conspicuous characteristic in a white-light image, an inter- 



568 planetary brightness has multiple origins such as a shock front and a CIR. 
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569 



These origins involve the compression of local plasma at the interface between 

570 two distinct streams. The shock could be an incident shock or a CME-driven 

571 shock. The CIR is formed as a result of the compression between the fast and 

572 slow streams, when both streams flow out of the rotating solar source surface 

573 at the same heliographic latitude. As a contrast, a shock is a transient distur- 

574 bance from solar eruptions, and a CIR is an ever-changing periodic structure 
in the background of interplanetary sp ace. Continu ously imaged in white 



575 



585 
586 
587 
588 



576 light, both a shock (IDavies et al.l . 120091 ) and a CIR ( IRouillard et al.l . 120081 1 



577 have the variance and movement in their optical brightness. Sometimes, the 

578 brightness of a CIR can be enhanced, when a preceding slow plasmoid is 

579 firstly swept and then entrained by the following fast CIR. Such a plasmoid 

580 imaged by an HI could be a pl asma blob disconnecte d from the cusp point 



58i of a coronal helmet st reamer (IRouillard et al.l . 120081 ) or a small-scale MC 



582 (IRouillard et al.l . 120091 ). Furthermore, when a CIR in the Sun-rooted spiral 

583 morphology blocks the trajectory of an energetic CME, the collision can lead 

584 to the CME becoming entrained by the CIR and the CIR being warped by 
the entrained CME. The interplanetary dynamics of the CME-CIR interac- 
tion would be manifested in white-light imaging as a more complex behavior 
of the brightness. Meanwhile, the IPS technique has its observational ca- 
pability to probe the micro-scale density fluctuation inside the macro-scale 

589 brightness imaged by an HI. Hence, the coordinated remote-sensing observa- 

590 tions of white-light and IPS are efficient to monitor the whole interplanetary 

591 space. 

592 The joint observations of white-light and IPS can provide the consistent 

593 observational evidence for the possible longitudinal deflection of a CME/shock 

594 in interplanetary space. The radial and latitudinal movements of a CME are 

595 recorded in the two-dimensional white-light image, once the CME is within 

596 the FOV of the HI and near the Thomson sphere surface of the HI. The lon- 

597 gitudinal movement of a CME could be inferred from the continuous white- 

598 light imaging with complementary IPS observation. The IPS signal gives the 

599 drifting speed of local density irregularity perpendicular to the IPS ray-path 
eoo from the Earth. Considered as the global bulk flow speed, the local IPS drift- 

601 ing speed derives the three-dimensional flow speed with the assumption of 

602 radial propagation. The derived radial flow speed can serve as a lower limit 

603 in the elongation of white-light imaging for an outwardly propagating CME, 

604 as interpreted in Section [5] and demonstrated in Figure [7b. The deviation 

605 of the predicted elongation-time curve suggests the non-radial propagation 

606 of the CME. If the latitudinal deflection is excluded from the HI imaging, 
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607 
608 



the non-radial propagation should come from the longitudinal direction. The 
longitudinal deflection can be again confirmed by the stereoscopic white-light 
609 imaging from the HI-A and HI-B instruments, as shown in Figure HJ For in- 
eio stance, if an Earth-directed CME is gradually deflected to the west and is 
en finally enclosed by the HI-A's Thomson sphere, the previously perfect sym- 

612 metry between the white-light images of two His is gradually broken, and 

613 the HI-A image becomes darker and darker. The deflection effect clarifies 

614 the disappearance of the CME-associated bright front in the FOV of HI-A. 

615 In fact, the shock aphelion in this paper does deviate to the west, because 

616 the shock front is quasi-perpendicular in the west and q uasi-parall el in the 



617 east as a result of the spiral configuration of the IMF (jHul . 11998 ). How- 

618 ever, the total deflection angle during the Sun-Earth space is only 3°, and 

619 is too small to be discerned by the observations of HI and IPS. The ignor- 

620 able longitudinal deflection in this paper is caused by the unimodal ambient 

621 stream of slow solar wind. If a CIR is incorporated into our model as one 

622 feature of the backgro und, the C ME deflection may be significant due to the 



623 CME-CIR interaction (jHul . Il998l ). Alternatively, if the initial eruptions of an 

624 early slow CME and a late fast CME are at an appropriate angular differ- 

625 ence, the contrary deflections of the two CMEs could be no t iceab le during 



626 the interplanetary process of oblique collision (IXiong et all 120091 ). These 

627 significant deflections result of the CME-CIR interacting or CME- 

628 CME coupling, and should be reflected from the observational signatures of 

629 white-light and IPS. These observational signatures will be further explored 

630 as a continuation to the preliminary results presented in this paper. 
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Figure 1: Density distribution n in the ecliptic plane for the propagation of an inter- 
planetary shock at (a) 19, (b) 23, and (c) 33 hours, following its initiation, with the 
interplanetary magnetic-field lines shown as solid black lines. The synthetic STEREO 
HI-A (HI-B) is at point 'A' ('B'), with its Thomson-scattering sphere marked by a white 
dotted circle, and field of view marked by two dash-dotted lines. Note that the synthetic 
STEREO A and B spacecraft are located at 45° in the longitude west and east from the 
Earth at point 'E', respectively. One solid line from the point 'E' denotes an IPS ray-path 
at the Earth with its elongation of 30°. Two solid lines from the point 'B' represent the 
elongation boundaries of relative-brightness enhancement (/ — Iq)/Iq > 0.4 in Figure [3] 
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Figure 2: Radial flow speed distribution v r in the ecliptic plane for the propagation of an 
interplanetary shock at (a) 19, (b) 23, and (c) 33 hours, following its initiation. 
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Figure 3: Relative-brightness enhancement (/ — Iq)/Io generated by the time-series as- 
sembling of an elongation slice within the ecliptic plane, according to synthetic data from 
the (a) HI-A and (b) HI-B. Here, Iq refers to the background brightness. 
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Figure 4: The density n (a) and speed v dli{t (b) profiles along an IPS ray-path for the 
background solar wind. Here, Udrift refers to the speed component perpendicular to the 
IPS ray-path. Below are (c) the Auto-Correlation Function (ACF) and (d) the Cross- 
Correlation Function (CCF) from the IPS observations. The baseline between the two 
radio antennas is 2000 km, aligned with the solar wind outflow. 
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Figure 5: The density n (a, e, i) and speed VdMt (b, f, j) profiles along an IPS ray-path 
for a shock propagation. Below are (c, g, k) the Auto-Correlation Function (ACF) and 
(d, h, 1) the Cross- Correlation Function (CCF) from the IPS observations. From the left 
to right, the three columns refer to three snapshots at 19, 23, and 33 hours, respectively, 
following the initiation of the event. 
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Figure 6: The profiles of the electron number n-r^ lB (a, d, g), radial bulk flow speed v r (b, 
e, h), and angular speed relative to the HI-B point ^hib (c, f, i) along a specific ray-path 
r H iB at the elongation of 26° at 19 hours, 30° at 23 hours, and 38° at 33 hours, viewed by 
the synthetic STEREO HI-B. 
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Figure 7: The time-dependent variances of (a) the radial aphelion r and (b) the propa- 
gation speed w S hock of a shock, as well as its manifestation in (c) the elongation of HI-B's 
field of view, shown as black solid lines. In panel (b), the bulk flow speed v r and the fast 
magnetosonic speed v r + Cf just downstream of the shock front are plotted as a dashed 
line and a dash-dotted line respectively. As a specific plasma parcel on the shock front is 
firstly detected as an IPS irregularity source at 19 hours, its later appearances in the radial 
distance from the Sun (a dashed line in panel a) and the elongation from the HI (a dashed 
line in panel c) could be estimated from the IPS observation at 19 hours together with 
two assumptions of radial direction and constant flow speed for the shock propagation. 
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Figure 8: The infrastructure of a multi-scale numerical model directly linking interplan- 
etary dynamics to corresponding observational signatures from white-light imaging and 
IPS. As a response to solar outputs, the macro-scale properties of electron number density 
n, bulk- flow speed V, magnetic field B, and temperature T in interplanetary space are 
firstly modelled from an MHD simulation. These macro-scale profiles along any ray-path 
can then generate synthetic white-light image and IPS data, using Thomson-scattering 
and radio-scattering theories respectively. 



31 



